ENERGY AND SYMMETRY OF ORDERED TWO-DIMENSIONAL DIPOLE 
CLUSTERS IN A PARABOLIC CONFINEMENT 
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We report an experimental study of two-dimensional dipole clusters formed by small magnetic 
particles floating at the air-water interface. The particles form hexagonally ordered clusters due 
to interplay between dipole-dipole repulsion and confining external forces. Optical images of the 
clusters allow determination of area S, total energy E, and the chemical potential /x for different 
particle numbers A'^ < 130. Asymptotic expressions for E{N), ^{N), S{N) account fairly well for the 
smooth part of our experimental results. However, on top of the smooth dependences we observe 
quasiperiodic fluctuations with dips at "magic" numbers corresponding to particularly symmetric 
cluster configurations. The magnitude and occurrence of these fiuctuations are related to the sym- 
metry of the cluster and to the center of mass position. This is in full agreement with the recent 
prediction of Koulakov and Shklovskii for Wigner crystal island. 

PACS numbers: 6f .46, 68.65+g, 36.40-c 

Confined two-dimensional (2D) systems of interacting 
particles, such as electrons in quantum dots vortices 
in plasma in superfluids Q and in mesoscopic su- 
perconductors adatoms; colloidal crystals at liquid 
interfaces electrons on the surface of the liquid he- 
lium 1^, have attracted much attention. These systems 
exhibit variety of interesting phenomena, including melt- 
ing P,[7|-[Tl|] , crystallization, shell structure ||l2|,|l3| and 
"magic" numbers p^ ]. While quantum mechanics is cer- 
tainly important here, many features of these systems 
arise solely from geometry and may be explained classi- 
cally. In particular, transport measurements on a quan- 
tum dot when it forms Wigner crystal at low electron con- 
centration or in a strong magnetic field, reveal periodic 
disappearance or merging of Coulomb blockade peaks 
. A classical explanation of this puzzling phenomenon 
has been proposed by Koulakov and Shklovskii jl^ who 
traced it to the intricate dependence of the crystal en- 
ergy on its symmetry. In particular, Ref. predicted 
that (i) the center of mass position in the unit cell of the 
2D Wigner crystal varies periodically with the number 
of particles; (ii) the center of mass sticks to positions of 
high symmetry. This prediction has been approved in 
computer simulations for hard-disc (short-range) inter- 
actions and Coulomb (long-range) potential so it is 
probably quite general and valid for all kinds of isotropic 
pair interactions. 

In this letter we demonstrate a macroscopic realization 
of 2D crystals, which allows to check experimentally the 
relation between crystal energy and symmetry. We use 
an array of permanent magnetic particles floating on the 
surface of the water |p6|-p8[. This array is confined by a 
nonuniform magnetic field of a circular coil with nearly 
parabolic distance dependence B = Bq{1 + / D^) [ p^ . 
Here, is the field in the center, and D is a param- 
eter (Fig.l). Similar systems consisting of superparam- 
agnetic colloidal particles on the liquid surface and in 



the presence of uniform magnetic field have been studied 
in the context of 2D-melting and pattern forma- 
tion Jiot . In contrast to these works, we use macroscopic 
particles for which the corresponding plasma parameter 
r ~ Eint/ksT K, 10^^* is enormous. Thus, the Brownian 
motion is insignificant. 

We study 2D cluster configurations vs particle number 
N . We add/remove one particle at the cluster bound- 
ary in the presence of the constant external field and 
wait for equilibration. To accelerate equilibration, we stir 
the particles using an oscillating magnetic field (Fig.l) 
with slowly decreasing amplitude (stimulated annealing) . 
After several minutes the particles self-assemble into a 
well-ordered cluster with quite reproducible configura- 
tion which remains stable for many days. 

An equilibrium particle configuration is a hexagonally 
packed circular cluster (Fig.l, lower panel) with pro- 
nounced shell structure. In all clusters there are six 
"topological" disclinations ||l^ that reside close to the 
boundary. This resolves the conflict between circularly 
symmetric magnetic confinement and hexagonal packing, 
arising from 1/r" interaction potential |20[| . In addi- 
tion, many clusters (and all clusters with N > 120) have 
few dislocations [|l0|,^ which reside at the boundary 
or are attached to the topological disclinations. There 
are highly symmetrical clusters with pronounced facets, 
which are usually free of dislocations. They occur at 
"magic" numbers (1,7,19,37,55,85. . . ) which roughly cor- 
respond to the number of particles in a perfect hexagon, 
namely, 3s(s -I- 1) -I- 1 where s = 0, 1, 2.... 

Cluster configuration is determined by the interaction 
energy Eint (dipole-dipole repulsion) and by the radially- 
dependent part of the confinement energy Efi^id, as fol- 
lows [H: 
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Here and rrii are the particle position and magnetic 
moment. For identical particles with magnetic moment 
m we introduce energy and distance scales: 

and recast the total energy in dimensionless form as: 

E = + Ef^^i, = ^ + ^ rf (3) 

Generally speaking, Eq.^ allows to find E in units of i?o 
(Eq.||) by summation of ^ N"^ terms using the measured 
particle positions. However, this may be reduced to the 
sum of only N terms due to special relation between Eint 
and Ef ieid in the equilibrium state. Indeed, let's consider 
an energy (Eq.|^) variation upon small uniform deforma- 
tion Ti — > ri(l + u), where u << 1. We find 

SE = -u{3E„,t - 2Ef,M) + u^{6E„,t + Ej.m) (4) 

In the equilibrium state the linear term vanishes which 
implies Eint ■ Efieid = 2:3. Therefore, 

E = ^Ef^eM = ll (5) 

Here, / = rf is the moment of inertia. It contains only 
N terms and is easily measured. Quadratic term in Eq.^ 
yields bulk compression modulus K = S^^&^E/dv? = 
6E/S, where S is cluster area. 

For each cluster with particle number N we take an 
image, using computer-controlled CCD camera, and de- 
termine the particle positions, using MATLAB. Then, we 
calculate the moment of inertia / = X^i''?; area S, and 
the center of mass position. The total energy E is calcu- 
lated through Eq.H with the accuracy of 0.15% as found 
in experiment with the same cluster which was shaked be- 
tween the measurements and then allowed to relax. The 
scale To (Eq.||) is found from the measurement of particle 
positions in the clusters consisting of 2-7 particles, whose 
configuration is evident. Figure 2 shows S{N), E{N), 
and the chemical potential fJ.{N) = dE/dN . These pa- 
rameters increase with N according to a power-law de- 
pendence ~ N^/^ which may be explained within a con- 
tinuum approximation. To this end we calculate total en- 
ergy in the limit of large N by replacing the sum (EqJ^) 
by the integral and assuming constant density p. We find 

where R is the cluster radius, 2a is the lattice constant, 
and p K. l/na^*. The first term in Eqj^ represents Ei„t, 
while the second term represents Efieu- Since according 
to Eq.| their ratio is 2:3, Eq.| yields = SR"^. Com- 
bining it with the obvious relation R w aN^, we find 



R — (1)^ A^io, a = (2^)^, and a single particle energy 
e = E/N = 4!^, so that 

e - 0.98A^*;^ - 1.57ivt;S' - 3.69iV^;if - 1.59A^ (7) 

Figure 2 shows that Eq.^ accounts fairly well for the 
experimentally found N^/^ dependence for 5, p, e (the 
prefactors are by 20% smaller and this will be discussed 
below) and for the linear dependence of K. Deviations at 
small N result from the surface energy and are accounted 
by a more advanced model (to be discussed elsewhere). 

Note deviations of S, p, e, K from smooth dependences 
(Fig. 2) which look like random fluctuations although 
they are quite reproducible and exceed the experimental 
errors by five times. These fluctuations appear in all clus- 
ter properties and are very prominent in p, where they 
amount up to 15%. Note that there are regions in which 
the chemical potential remains flat or even decreases with 
addition of new particles (Fig. 2). 

Figure 3 shows expanded view of the fluctuations in 
total energy SE (after subtraction of a smooth back- 
ground). It also shows second derivative of the total 
energy, g^^f* . In the context of quantum dots the lat- 
ter corresponds to charging energy. Pronounced positive 
peaks in ^^^f occur for "magic" clusters which have al- 
most perfect symmetry. The energy of these clusters is 
lower (note dips in 6E at Fig. 3) while the density and 
the bulk compression modulus is higher (not shown here). 
The magnitude of the fluctuations is 6E « O.le and 
SS ~ S/N. In three experimental runs (one for decreas- 
ing N and two for increasing N) the energy dips occur 
almost at the same magic numbers (7,19,35,55,86,118..). 
While magic numbers are well-known for free-standing 
3D clusters [ p^ , and almost the same numbers have been 
found in computer simulations for 2D clusters ||l^,|l^,|l5| , 
we are unaware of previous experimental observations of 
magic 2D clusters in lateral conflning fleld. 

Koulakov and Shklovskii ||l^ have recently related 
magic numbers and energy fluctuations to the cluster 
symmetry and to the center of mass position in the crys- 
talline lattice. While the absolute position of the center 
of mass is always at the minimum of the external po- 
tential, its position with respect to the unit cell is not 
flxed. Figure 4 shows a unit cell of the hexagonal lattice 
where (A), (B) and (C) represent points of high symme- 
try (six-fold, three- fold and two- fold, respectively). Ref. 
|]l5| shows that the energy minimum is achieved when the 
center of mass of the whole cluster is located at one of 
these points. We measured distribution of the center of 
mass positions in the unit cell of our 2D dipolc clusters 
and indeed found that it is strongly peaked at (A), (B) 
or (C) (Fig. 4) . A lower panel in Fig. 3 also shows Ncm, a 
number of nearest neighbors at the center of mass posi- 
tion vs particle number. [We take one of the points (A), 
(B), or (C) which is the closest to the center of mass and 
assign Ncm = 6 for (A), Ncm = 3 for (B), and Ncm = 2 
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for (C)]. Ncm changes quasiperiodically in N^. There 
are continuous ranges where the center of mass resides 
in the position of six- fold symmetry (A). In these ranges 
the dips of energy and pronounced peaks of occur. 
Shallow dips correspond to position (B). Ref. i|5| shows 
that upon addition of new particles, the cluster undergoes 
energetically expensive elastic and plastic (formation of 
dislocations) deformations in order to keep the center of 
mass at the position of high symmetry. This occurs until 
another position of high symmetry becomes energetically 
favorable. Then, upon addition of one more particle, an 
avalanche occurs. This means rearrangement of many 
particles, whereby the center of mass moves to a large 
distance a. The avalanches take place with periodicity 
of ~ iVa which corresponds to addition of a new crys- 
talline row. In Fig. 4 we present the distance between the 
center of mass and the nearest particle and indeed ob- 
serve that it does not vary gradually, but in a step-like 
fashion. 

We estimate the magnitude of fluctuations as fol- 
lows. Assume a perfectly symmetric cluster. Addition 
of new particles does not lead to the center of mass 
displacement until approximately one crystalline row or 
A7V ~ (Af/3)i/2 particles are added. If the cluster 
were incompressible, corresponding area increase would 
be SS ~ SAN/N. Since our clusters are compressible, 
immobility of the center of mass means that the cluster 
adopts new particles without appreciable change in area. 
Therefore, deviation from the smooth dependence S{N) 
is SS ~ S/{3Ny/'^ and it arises from cluster compress- 
ibility and deformation. If the deformation were purely 
elastic, then u — SS/2S and the corresponding energy 
change would be 6E = KSv? /2 ^ 0.25e. Experimentally 
observed energy fluctuations (O.le— 0.05e) are smaller due 
to formation of the dislocations that lead to relaxation 
of elastic stresses. 

Our experimental results should be compared to nu- 
merical simulations of Ref. ||l^ for 2D dipole clusters 
{N < 80) in a parabolic confinement. Experimentally 
found shell structures for < 37 are almost identical to 
those found in simulations, while at 37 < A^ < 80 they 
differ in half of the cases (the difference is mostly in the 
configuration of outer shells) and are more close to those 
found in numerical simulations for Coulomb clusters p2[ . 
Energy fluctuations found in experiment and in Ref. 
occur almost at the same "magic" numbers, although 
the experimentally found magnitude of SE is four times 
bigger than that found in simulations |jl^. The single 
particle energies at small A^ almost do not differ, while 
at large A" the experimental values are lower than those 
found in numerical simulations (see Fig. 3). We find em- 
pirical dependence, Sexp/^sim 1 — 0.1 log A^ in the whole 
range of A^. This small difference between eexp and egim 
may be attributed to (i) additional attractive capillary 
forces between the particles which logarithmically decay 



with distance |T^; (ii) deviation of the confining field 
from parabolicity which appears mostly at large A^ when 
the cluster radius becomes comparable to the radius of 
the coil. Both these factors lead to the cluster compres- 
sion and to the underestimate of the total energy using 
Eq.|. 
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FIGURE CAPTIONS 



Fig.l Experimental setup. Small rare-earth permanent 
magnets encapsulated within styrofoam discs are 
floating on the surface of the water and repel each 
other. They are confined within a nonuniform mag- 
netic field produced by external coil. Container di- 
ameter is 100 cm, magnet size is 3 mm. Lower panel 
shows two cluster images with superimposed Delau- 
nay triangulation. Gray circles stand for normal 
coordination number (Z = 6), filled circles stand 
for Z = b (positive disclination), and open circles 
stand for Z = 7 (negative disclination). A^=85 - 
a "magic" cluster. Note hexagonal structure and six 
topological disclinations at the boundary. A^=70- 
an ordinary cluster. Note a dislocation (5-7 pair) 
attached to a disclination. 

Fig. 2 Single particle energy e, the chemical potential fi 
and cluster area S vs particle number N (in dimen- 
sionless units). Continuous lines show experimen- 
tal data (smoothed); dashed lines show asymptotic 
A^'^/^dependence; dashed-dotted line shows results 
of numerical simulations of Ref. ||l3| (multiplied by 
0.8) for the energy of 2D dipole clusters in parabolic 
confinement. The inset shows bulk compression 
modulus K. 

Fig. 3 Fluctuations of the total cluster energy 6E and 
"charging energy" |^ (note inverted scale) vs par- 
ticle number N. The lowest curve shows number 
of nearest neighbors, Ncm, to the center of mass 
position. Note correlation between three curves. 
Big positive values of g^^f (dips in the figure) 
correspond to dips in SE. Deep dips occur when 
Ncm = 6 and shallow dips correspond to Ncm — 3. 
"Magic" numbers indicated on the Figure corre- 
spond to energy dips and almost perfect clusters. 

Fig. 4 Center of mass position, (a) Unit cell and the posi- 
tions of high symmetry, (A),(B), and (C). (b) Dis- 
tribution of the center of mass position in the unit 
cell. Each point corresponds to a cluster with a cer- 
tain particle number N. Note sticking of the center 
of mass to positions (A),(B), and (C). Lower panel 
shows the distance between the center of mass and 
the nearest particle. Note periodicity. 
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This figure "fig2.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0008 117vl 



This figure "fig3.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0008 117vl 



This figure "fig4.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0008 117vl 



This figure "figl.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0008 117vl 



